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Abstract 



Coulomb interaction between electrons on p-orbitals of oxygen atom in 



^j>. I strongly correlated compounds is not negligible, since its value (Up) has com- 

^^ I parable order of magnitude with the value of Coulomb interaction on d- 



orbitals of transition metal atom (Ud)- We investigate the effect of taking 



o 
o . 

C^ , into account Coulomb correlations in oxygen p-shell in addition to the corre- 



ct ] lations in the transition metal d-shell in frame of the LDA+U method. Our 

, , calculations for NiO, MnO and La2Cu04 show that this additional correction 

Q I in general improves the agreement with experimental data for the spectral 

O ' 

(energy gap values, relative position of the main peaks in X-ray photoemis- 

K> I sion spectroscopy (XPS) and Bremsstrahlung isohromate spectroscopy (BIS)) 

and magnetic properties (magnetic moment values and intersite exchange in- 
teraction parameters values). 

71.27.+a Strongly correlated electron systems; heavy fermions - 74.25.Jb Elec- 
tronic structure - 79.60.-i Photoemission and photoelectron spectra 
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I. INTRODUCTION 

During last several decades the electronic structure calculations from first principles 
became an important part of the solid state theory. The solution of such essentially many- 
particle problem as calculation of band structure of real materials is impossible without 
rather severe approximations. The most famous and commonly used approximation in 
ab initio electronic structure calculations is density functional theory (DFT) |I| within the 
local spin density approximation (LSDA). But, as an approximation based on the homo- 
geneous electron gas theory ^, the LSDA is valid only for compounds with slow varying 
through the crystal charge density. In other words, the LSDA must describe well only 
the delocalized electronic states (broad bands). Nevertheless the LSDA is sometimes able 
to give correct ground state properties for the systems with rather narrow bands (see for 
example 0). 

The most unusual physical properties were found in the systems with strong electron- 
electron Coulomb correlations (such as Mott insulators, high-Tc superconductors, etc.). 
These systems have been intensively investigated during last 20 years both by experimen- 
talist and theoretician communities. All intriguing features of these systems come from 
the existence of nearly localized electronic states (narrow bands) such a.s d or f states of 
transition metals ions or rare-earth metals ions, respectively. To mark out the localized 
state one can apply the following criteria: kinetic energy of localized states is of the same 
or even smaller order of magnitude as the energy of Coulomb interactions. It is known that 
for strongly-correlated systems the LSDA often fails (high-Tc related compound La2Cu04; 
insulating, antiferromagnetic transition metal oxides). However constrained LSDA calcula- 
tions 0] give Coulomb interaction parameter values in surprisingly good agreement with the 
experimental estimations P- |ll|] . 

Several approaches were built on the LSDA basis repairing its deficiency in describing 
Coulomb interaction between localized states. The most popular methods are the self- 



interaction correction method (SIC) [0,0 and the LDA-I-U method [0. 



The basic problem of the LSDA is the orbital-independent potential which does not al- 
low to reproduce Coulomb interaction derived energy splitting between occupied and empty 
subbands. The SIC method solves this problem by introducing the orbital-dependent po- 
tential correction which explicitly substracts the self-interaction present in the LSDA. This 
method restores correct electronic properties of the transition metal oxides where the LSDA 
fails. However the self-interaction correction for the d states is so strong, that when one im- 
plements SIC potential only to the d-orbitals of transition metal then the oxygen p-orbitals 
do not shift from the LSDA obtained positions and the occupied d-band lies much lower in 
energy than oxygen valence band, which does not agree with the spectroscopy data. How- 
ever the values of energy gaps and the spin magnetic moments are in rather good agreement 



with experiment ||T3[. To improve this situation one can treat all valence states (namely the 



transition metal ions (i-orbitals and oxygen ^-orbitals) as localized and apply SIC potential 



to all of them [|lg]. In this case structure of occupied bands is well reproduced, but the value 



of energy gaps will be overestimated |T^ 



Another way to overcome the well known disadvantages of the LSDA is the LDA+U 
method, which gives better agreement with experimental spectra |]I1| . The LDA+U method 



corresponds to the static limit of recently developed new many-body approach — the dy- 
namical mean- field theory (DMFT) [|l^. In its standard form the LDA+U takes into account 
only Coulomb interaction between d (or /) electrons of transition metal ions. In the present 
paper we investigate the problem of Coulomb interaction between oxygen p electrons and 
show that the inclusion of the corresponding term in LDA+U equations leads to signifi- 
cant improvement of agreement between calculated and experimental spectral and magnetic 
properties. 

II. METHOD OF CALCULATION 

The main idea of the LDA+U method is to add to the LSDA functional the term E^ 
corresponding to the mean-field approximation of the Coulomb interaction in multiband 



Hubbard model. 

where p'^(r) is the charge density for spin-cr electrons and E^^^^\p'^{j)\ is the standard 
LSDA (Local Spin-Density Approximation) functional. Eq. (|l|) asserts that the LSDA is 
sufficient in the absence of orbital polarizations, while the latter are driven by, 

E^{{n''}\ = lE{m},a{{m,m" \ Vee I m',m'")n^^,n;;,^^,„ + 

(2) 

[{m^m" I Vee I m' ,m"') — {m,m" \ Vee \ f^"' if^'))'Knm'''^'rn"m"'}i 

where Vee are the screened Coulomb interactions among the d electrons. Finally, the last 
term in Eq. (|I]) corrects for double counting (in the absence of orbital polarizations, Eq. (|1]) 
should reduce to E^^^^) and is given by 

Edc[{n'']] = hjN{N - 1) - ]^J[N\N^ - 1) + N^N^ - 1)], (3) 

were N'^ = Tr(n'^.^,) and A^ = A^^ + A^-'-. U and J are screened Coulomb and exchange 
parameters [p|,p!0|. 

In addition to the usual LSDA potential, an effective single-particle potential to be used 
in the effective single-particle Hamiltonian has the form: 

H = Hlsda + Yl I inlma)V^^,{inlm'a \ (4) 

mm' 
^mm' — Sm"m"'{('^5'"^ | Ke | ^ 5 ""^ )^rn/'m"'~^ 

{{m,m" I Vee I m',m"') — {m,m" \ Vee \ m'" ,m'))n'^„^Ni}— (5) 

U{N -\) + J{N'' -\). 

The matrix elements of Coulomb interaction can be expressed in terms of complex spherical 
harmonics and effective Slater integrals F^ |jT^ as 



(m, m" I Vee I m\ m'") = ^ ak{m, m', m", m"')F^, (6) 

k 

where < A; < 2/ and 

47r '^ 
au{m,m\m'\m"') = —— ^ (Zm | n, | /m')(/m" | F,; | Zm'") 

For d electrons one needs F°, F"^ and F^ and these can be hnked to the Coulomb- and 
Stoner parameters U and J obtained from the LSDA-supercell procedures via U = F^ 
and J = [F"^ + F^)/1A, while the ratio F'^/F^ is to a good accuracy a constant ~ 0.625 
for the 3d elements [^,^]. (For / electrons the corresponding expression is J = (286F^ + 
195F^ + 250F^)/6435). The Coulomb parameter U is calculated as a second derivative of 
the total energy (or the first derivative of the corresponding eigenvalue) in respect to the 
occupancy of localized orbitals of the central atom in a supercell with fixed occupancies on 
all other atoms [§]. 

If one neglects the exchange and non-sphericity of the Coulomb interaction (which is 
exact in the case of the fully occupied or empty band) the potential correction will have the 
more simple form: 

V^ = U{^ - nO (7) 

where rij is the occupancy of z-orbital. Then for fully occupied state LDA+U potential 
correction would be the shift to the lower energies on U/2, while for empty states it gives 
an upward shift on the same value. So the LDA+U gives correct splitting between occupied 
and empty subbands equal to the Coulomb interaction parameter U. 

In the LDA+U approach the Coulomb interactions are taken into account conventionally 
only on d-orbitals of transition metals. However it is known that Coulomb interactions 
between electrons on ;?-orbitals of oxygen have comparable order of magnitude 1^,0 with 
the corresponding d — d Coulomb interactions and so must be taken into consideration on the 
same footing as for c?-orbitals. The usual justification for omitting of U on oxygen ^-shell is 
that the oxygen shell is fully occupied and the correlation effects between electrons (or rather 



holes) in it can be neglected due to the small number of holes in ground state. However the 
LDA+U equations (|^) will give nonzero correction for the fully occupied oxygen band: 

Vp = -Up/2 (8) 

This potential correction must be applied to the orbitals forming oxygen band, however 
corresponding Wannier functions (in contrast to d states) are far from being of pure 0{2p) 
character because they have very strong admixture of s and p states of transition metal 
ions and other extended orbitals. Since the main influence on the electronic structure is 
the change of the energy separation between the oxygen p-band and the transition metal 
d-band, the upward shift in energy of the transition metal d-band on Up/2 will be equivalent 
to the shifting down of the oxygen j)-band on the same value. Thus in our calculations we 
added Up/ 2 term to the diagonal matrix elements of the LDA+U potential correction (^. 
We call this extension of the LDA+U method in the paper as the LDA+U^'^^^-' . 

Recently the modified LDA+U*^"'"''^-' method was used by Korotin et al. for investigation 
of charge and orbital ordering effects in Lay/gSri/gMnOa compound |2^. The inclusion of 



Coulomb interactions in the oxygen j)-shell was found to be crucial in that calculation, since 
it controls the value of charge transfer energy between Mn(3(i) and 0{2p) valence states and 
significantly enhances the tendency of localization in this system. 

In the present paper we report the results with the use of the modified LDA+U^'^^'^-' 
method for the typical strongly correlated transition metal oxides NiO, MnO and La2Cu04. 
We show that inclusion of the correlations in oxygen j>-shell leads to the better agreement 
with the experimental data for the main peaks position in X-ray photoemission spectroscopy 
(XPS) and Bremsstrahlung isohromate spectroscopy (BIS) spectra in comparison with con- 
ventional LDA+U calculated spectra. Not only spectral properties, but both spin magnetic 
moments and intersite exchange interaction parameters J^x for NiO, MnO and La2Cu04 are 
in better agreement with the corresponding experimental data. 



III. RESULTS AND DISCUSSION 

The important part of the LDA+U calculation scheme is the determination of Coulomb 
interaction parameters U and J in equations (^): Coulomb parameter Up for ^-orbitals of 
oxygen, Ud for transition metals ion and Hund's parameter J for (i-orbitals of transition 
metals. To get Ud and J one can use the supercell procedure [rU,rT| or the constrained 



LSDA method [Q, which are based on calculation of the variation of the total energy as a 
function of the local occupation of the (i-shell. We took the values of Ud and J parameters 
listed in Table I previously calculated in WM . The problem is how to determine the Coulomb 



parameter Up. 

Due to more extended nature of the 0{2p) Wannier states in comparison with transition 
metal d states, the constrained occupation calculations can not be implemented as easy as for 
the d-shell of transition metals. Nevertheless several independent and different techniques 
were used for this purpose previously by different authors. McMahan et al. estimated the 
value of Up in high-Tc related compound La2Cu04 using the constrained LDA calculation 
where only atomic-like 0(2p)-orbitals within oxygen atomic spheres were considered instead 
of the more extended Wannier functions. The corresponding value of Coulomb interaction 
parameter Up was obtained as 7.3 eV . This value can be considered as the upper limit of 
the exact Up. The LDA calculations gave the estimation that only 75% of Wannier function 
density lies in the oxygen atomic sphere so that renormalized value of Coulomb interaction 
parameter for oxygen Wannier functions is Up = (7.3) x (0.75)^ = 4.1 eV 0. 

Later Hybertsen et al. suggested the scheme to calculate Up, which consists of two steps: 
(i) via constrained-density-functional approach one can obtain the energy surface E{Nd, Np) 
as a function of local charge states and (ii) simultaneously extended Hubbard model was 
solved in mean-field approximation as a function of local charge states Nd and Np. Cor- 
responding Coulomb interaction parameters were extracted as those which give the energy 
surface matching the microscopic density- functional calculations results 0. The obtained 
values for Up are 3 -i- 8 eV depending on the parameters of calculations. 



Another way to estimate Up is to use Auger spectroscopy data, where two holes in 0(2p)- 
shell are created in the excitation process. Such fitting to the experimental spectra gave the 
value of Up = 5.9 eV [gSf. In our LDA+U^'^+p) calculations we used Up = 6 eV. 

To solve the LDA+U hamiltonian we implemented the self-consistent tight- 
binding (TB) linear muffin-tin orbitals method (LMTO) in the atomic sphere approxima- 
tion (ASA) p^- p^] . For the calculations we choose the classical strongly correlated transition 
metal oxides NiO, MnO and La2Cu04, which were well investigated by experimental and 
theoretical methods. 

Comparison between the LDA+U (left column) and the LDA+U'-^^*'-* (right column) cal- 
culated density of states (DOS) of NiO, MnO and La2Cu04 is presented in figures 1, 2 and 3. 
For all compounds one can see that the main difference between the LDA+U*-'^'^^'' and the 
LDA+U calculated densities of states is the increased energy separation between the oxy- 
gen 2p and transition metal 3d bands. The larger value of "charge transfer" energy (0(2p)- 
Me(3d)) (Me=Ni,Mn,Cu) leads to the enhanced ionicity and decreased covalency nature of 
the electronic structure: the unoccupied bands have more pronounced 3d character and the 
admixture of oxygen states to those bands becomes weaker. 

The ground state is correctly described both by LDA+U and LDA+U*-'^"'"^'' calculations 
as antiferromagnetic insulator for all compounds. The values of energy gaps |^ and spin 



magnetic moments are presented in tables II and III. One can see that the values obtained 
in the LDA+U*^'^^^^ calculations are in general in better agreement with experiment than the 
LDA+U calculated values. While the increasing of the energy gap values with applying Up 
correction was obviously expected with the increasing of "charge transfer" energy in the 
compounds belonging to the class of " charge transfer" insulators p^ , the increasing of the 



magnetic moments values is more complicated self-consistency effect due to the increased 
ionicity in the LDA+U*-'^"*"^-' calculations comparing with the LDA+U results. 

In Fig. 4 the DOS obtained by LDA+U^'^+p^ and LDA+U calculations for MnO and NiO 
compounds are compared with the superimposed XPS and BIS spectra corresponding to the 
removal of an electron (the occupied bands) and addition of an electron (the empty bands), 

8 



respectively. The better agreement with the experimental data of position of the main peaks 
of unoccupied band relative to the occupied one is the direct confirmation of the importance 
of taking into account Coulomb interactions in oxygen 2j)-shell. 

The values of the intersite exchange interaction parameters Jex depend on the param- 
eters of the electronic structure in a rather indirect implicit way. The developing of the 
good calculating scheme for exchange parameters is very important because the ab-initio 
calculation is often the only way to describe the magnetic properties of complicated com- 
pounds such as for example "spin-gap" systems |^. Recently Solovyev et al. |^ did very 
through analysis of the exchange interaction parameters for MnO calculated using different 
methods of electronic structure calculations. They used the positions of the Mn(3(i)-spin-up 
and Mn(3(i)-spin-down bands relative to the oxygen 2p states as adjustable parameters to 
fit the values of exchange interaction for the nearest and second Mn-Mn neighbors. Their 
results gave nearly the same splitting between Mn(3(i)-spin-up and Mn(3(i)-spin-down states 
as in standard LDA-I-U calculations (10.6 eV) but the position of those states relative to the 
oxygen band was shifted approximately on 3 eV up relative to the LDA+U case. It is practi- 
cally the same as we have in our LDA+U*^'^^^'' calculations, because with Up = 6 eV the shift 
of the position Me(3(i)-band relative to the oxygen 0(2p)-band is equal to Up/ 2 = 3 eV. 

Comparison between LDA-I-U and LDA-I-U*-'^^^'' calculated Jex parameters and experi- 
mental data is presented in table IV. Jex were calculated from Greens function method as 
second derivatives of the ground state energy with respect to the magnetic moment rotation 



angle |]31| , |32[| . Again one can see that in general the LDA-I-U*^'^'^^^ gives better results than 



the LDA-I-U, especially for MnO compound. 

IV. CONCLUSION 

The method for inclusion of Coulomb interactions between oxygen p electrons in the 
calculation scheme of the LDA-I-U method was proposed. The main effect was found to be 
the increasing of "charge transfer" energy parameter (the separation of 0(2^) and Me(3t?) 



states). As the result, the spectral and magnetic properties of the typical strongly cor- 
related transition metal oxides NiO, MnO and La2Cu04 were found in better agreement 
with experimental data than in the conventional LDA+U method where only correlations 
between Me{3d) state are taken into account. 
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FIG. 1. La2Cu04 DOS calculated by the LDA+U (left column) and the LDA+U^'^+p) (right 

column) methods. (On all figures the total DOS is presented per formula unit, the DOS of particular 

states are per atom. Fermi energy corresponds to zero.) 
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FIG. 2. MnO DOS calculated by the LDA+U (left column) and the LDA+U^'^+p) (right column) 
methods. 
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FIG. 3. NiO DOS calculated by the LDA+U (left column) and the LDA+U^'^+p) (right column) 

methods. 
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FIG. 4. DOS calculated by the LDA+U (dashed line) and the LDA+U^'^+p) (sohd line) Ni(3(i) 
and Mn(3d) in comparison with superimposed XPS and BIS spectra. 
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Spectroscopic data: 
MnO: G.A. Sawatzky and J.W. Allen, Pliys. Rev. Lett. 53, 2339 (1984). 
NiO: J. van Elp, R. H. Potze, H. Eskes, R. Berger, and G. A. Sawatzky, Phys.Rev.B 44, 
1530 (1991). 
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TABLE I. Coulomb parameters 


Ud 


TABLES 

and Hund's parameters 


J [eV) used in calculations. 


La2Cu04 




MnO 


NiO 


Ud 8.0 
J 1.0 




6.9 
0.86 


8.0 
0.95 


TABLE II. Calculated and experimental values of 


energy gaps {eV). 


LDA+U 




LDA+U^'^+P) 


Experiment 


La2Cu04 0.7 
MnO 3.8 
NiO 1.8 




2.0 
4.5 
2.8 


2.0^^ 

3.6-3.8^ 

4.3^4.0'=' 
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TABLE III. Calculated and experimental values of spin magnetic moments (pb)- 





LDA+U 


LDA+U^'^+P) 


Experiment 


La2Cu04 


0.45 


0.68 


0.60'^ 


MnO 


4.51 


4.59 


4.79^4.58^ 


NiO 


1.50 


1.64 


1.77*^, 1.64^1.90-^ 



''Y. Endoh, et ai, Phys. Rev. B 37, 7443 (1988). 
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TABLE IV. Calculated and experimental values of intersite exchange interaction parame- 
ters Jex (meV). 



LDA+U 



LDA+U^'^+P) 



Experiment 



iLa2Cu04 


82.9 


100.9 


136 ± 5" 


2MnO 


5.4 9.3 


5.4 5.1 


4.8^5.4^ 5.6^5.9'^ 


2NiO 


0.8 23.2 


0.2 19.4 


1.4^^ 19.0'=' 



^Cu-Cu exchange parameter between nearest Cu atoms in plane. 

^Me-Me exchange parameters between nearest and second neighbors. 
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